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SUMMARY 

Radioactivity  induced  in  detectors  by  protons  and  secondary  neutrons 
limits  the  sensitivity  of  spaceborne  ganma-ray  spectrometers.  Three  dimensional 
Monte  Carlo  transport  codes  have  been  employed  to  simulate  particle  transport  of 
cosmic  rays  and  inner-belt  protons  in  various  representations  of  the  Gamma  Ray 
Observatory  Spacecraft  and  the  Oriented  Scintillation  Spectrometer  Experiment. 
Results  are  used  to  accurately  quantify  the  contributions  to  the  radioactive  back¬ 
ground,  assess  shielding  options  and  examine  the  effect  of  detector  and  space¬ 
craft  orientation  in  anisotropic  trapped  proton  fluxes. 
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PARTICLE  TRANSPORT  SIMULATION  FOR  SPACEBORNE,  Nal  GAMMA-RAY  SPECTROMETERS 

*  Dr  C  S  Dyer,  P  R  Truacott,  A  J  Sims 
tDr  C  Comber,  Dr  N  D  A  Hammond 


Abstract 

Radioactivity  induced  in  detectors  by  protons  and  secondary  neutrons  limits 
the  sensitivity  of  spacebome  gamma-ray  spectrometers.  Three  dimensional  Monte 
Carlo  transport  codes  have  been  employed  to  simulate  particle  transport  of 
cosmic  rays  and  inner-belt  protons  in  various  representations  of  the  Gamma  Ray 
Observatory  Spacecraft  and  the  Oriented  Scintillation  Spectrometer  Experiment. 
Results  are  used  to  accurately  quantify  the  contributions  to  the  radioactive 
background,  assess  shielding  options  and  examine  the  effect  of  detector  and 
spacecraft  orientation  in  anisotropic  trapped  proton  fluxes. 

Introduction 

Radioactivity  induced  in  scintillation  detectors  by  cosmic  rays,  trapped  protons 
and  secondary  neutrons  produces  an  important  source  of  background  in  spacebome 
gamma-ray  spectrometers*.  For  the  early  detectors,  which  employed  relatively 
lightweight  crystals  of  Nal  and  Csl  (a  few  pounds)  on  modestly  sized  spacecraft 
(a  few  hundred  pounds),  typical  particle  pathlengths  through  the  spacecraft  and 
detector  material  were  significantly  less  them  the  nuclear  interaction  length 
and  the  direct  spallation  interactions  of  primary  particles  sufficed  to  explain 
the  observed  radioactivity*.  For  the  next  generation  of  instruments  to  be 
carried  on  the  Gamma  Ray  Observatory  individual  detector  masses  of  order  2  tons 
are  to  be  carried  on  a  15  ton  spacecraft.  Typical  particle  pathlengths  exceed 
the  nuclear  interaction  length  and  the  consequent  importance  of  secondary 
particles  requires  the  application  of  particle  transport  codes.  In  a  previous 
paper 2  1-D  transport  calculations ,  which  employed  the  ANISN  *  neutron  code  in 
spherical  geometry,  were  used  to  assess  the  importance  of  secondary  neutron 
capture  on  Nal  and  explore  the  efficacy  of  ‘LiH  shield  layers  around  the  central 
Nal  elements.  Such  layers  were  in  fact  shown  to  be  counterproductive  and  this 
result  was  confirmed  by  preliminary  results  from  a  3-D  simulation  of  trapped 
proton  transport  in  a  more  representative  geometry*.  In  the  current  paper  the 
full  3-D  simulation  is  applied  to  a  number  of  detector  and  spacecraft  configur¬ 
ations  in  order  to  assess  the  dependence  of  the  various  background  components  on 
spacecraft  and  detector  geometries  and  particle  anisotropies. 

Components  of  the  Calculation 


The  Transport  Suite 


The  codes  used  in  this  study  all  employ  Monte  Carlo  techniques  and  include  the 
High  Energy  Transport  Code  (HETC)S  for  protons  and  energetic  hadrons  and  the 
MORSE*  code  for  neutrons  of  less  than  15  MeV. HETC  is  a  nucleon/meson  transport 
code  which  uses  an  intranuclear  cascade  and  evaporation  model7  to  simulate  the 


*  Space  Department,  Royal  Aerospace  Establishment,  Farnborough,  Hampshire,  England, 
t  SD-Scicon  Pic,  Abbey  House,  Farnborough  Road,  Farnborough,  Hampshire,  England. 
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nuclear  interactions.  MORSE  is  a  coupled  neutron/gamma  ray  code  which 
utilises  multigroup  scattering  and  absorption  cross-section  data.  It  is 
used  here  to  transport  neutrons  over  the  energy  range  from  15  MeV  down  to 
thermal  energies  using  the  neutron  source  generated  by  HETC .  These  codes 
share  a  common  combinatorial  geometry  package.  Other  codes  used  to  predict 
spallation  product  yields,  time  variations  in  induced  radioactivity  and  energy- 
loss  spectra  for  decay  products  are  as  discussed  in  Ref  1. 

Detector  and  Spacecraft  Geometry 

The  detectors  and  spacecraft  under  study  are  representations  of  the  Oriented 
Scintillation  Spectrometer  Experiment  (OSSE)*  to  be  carried  on  the  Gamma  Ray 
Observatory.  The  OSSE  detector  employs  four  identical  telescopes,  each  of 
which  comprise  a  13  inch  diameter  central  Nal/Csl  phoswich  (a  phosphor  sandwich 
with  a  4  inch  thickness  of  Nal  plus  3  inches  of  Csl  and  utilising  the  different 
pulse  shapes  in  the  two  materials  to  veto  events  arriving  from  the  rear) , 
provided  with  an  active  Nal  shield  of  about  3  inch  thickness  and  a  passive 
slat  collimator  of  tungsten  to  give  an  opening  angle  of  5°  x  11°.  The 
telescopes  are  orientable  through  about  192°  in  pairs.  The  combinatorial 
geometry  representation  is  illustrated  in  Figure  1.  For  simplicity  the  tungsten 
collimator  is  represented  by  a  diffuse  mass  of  lower  density  having  the  same 
total  mass  as  the  actual  slat  arrangement.  The  ‘LiH  layer  around  the  inner 
phoswich  was  given  a  thickness  of  0.3  cm  to  represent  the  ‘LiF  impregnated 
silicone  actually  incorporated  in  the  instrument.  This  material  was  incorporated 
in  order  to  reduce  the  influence  of  thermal  neutrons  arising  from  hydrazine 
tanks  and  plastic  scintillators  carried  on  the  spacecraft.  In  certain  simulations 
this  thickness  has  been  increased  up  to  5  cm  to  investigate  the  efficacy  in 
reducing  the  more  energetic  neutrons  from  the  spacecraft  and  outer  Nal  shield. 

Such  increased  thicknesses  proved  to  be  counterproductive  in  this  case*.  For 
certain  simulations  adjacent  spacecraft  material  was  represented  by  10  cm  of 
aluminium  on  two  sides  as  in  Figure  1.  In  other  cases  an  isolated  telescope 
was  studied  and  in  further  cases  a  diffuse  aluminium  model  of  the  spacecraft 
was  employed  (Figure  2).  Intercomparison  of  these  cases  enables  consideration 
of  the  relative  importance  of  secondaries  arising  from  the  detector,  local 
material  and  spacecraft. 

Proton  Fluxes 


The  incident  trapped  proton  and  cosmic  ray  fluxes  and  their  spectra  are  those 
predicted  for  a  28.5°,  450  km  orbit  in  1991.  As  discussed  in  Ref  4,  trapped 
protons  have  a  characteristic  energy  of  100  MeV  and  probably  do  not  exceed 
600  MeV,  while  the  majority  of  cosmic  ray  interactions  arise  from  protons  in 
the  energy  range  from  6  to  30  GeV.  The  exact  flux  of  trapped  protons 
experienced  in  the  South  Atlantic  Anomaly  remains  in  some  doubt,  in  part  due 
to  uncertainties  in  the  models  and  in  part  due  to  debate  as  to  the  correct 
procedures  for  incorporating  the  influence  of  the  time  varying  geomagnetic 
field.  In  the  light  of  this  debate  Stassinopoulos’  has  recently  revised  the 
GRO  proton  fluence  estimates  and  uses  two  alternative  approaches  to  arrive  at 
integral  daily  fluences  in  excess  of  100  MeV  of  3.6  x  10s  and  7.1  x  10*  cm-1  . 

In  this  paper  we  therefore  use  a  value  of  5  x  10*  with  an  uncertainty  factor  of 
two.  This  fluence  is  a  factor  of  three  lower  than  the  value  used  in  Ref  4. 

For  simplicity  isotropic  incidence  has  been  used  in  most  cases.  However  use 
of  these  3-D  codes  enables  the  influence  of  alternative  distributions  to  be 
studied.  Such  effects  are  likely  to  be  more  important  for  trapped  protons  as 
they  are  highly  anistropic  and  also  less  penetrating  than  cosmic  rays.  Recent 
studies  of  Space  Station  shielding  have  led  to  reconsideration  of  the  spatial 
distribution  of  South  Atlantic  Anomaly  protons1 *  .  Using  this  work  proton 
distributions  have  been  generated  as  illustrated  in  Figure  3,  which  shows  that, 
in  the  SAA,  protons  are  close  to  their  mirror  point  and  have  pitch  angles  close 
to  90°.  There  is  a  further  East-West  anisotropy  which  increases  with  energy  and 
is  due  to  the  removal  of  particles  by  the  atmosphere. 
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Results 


Sample  results  obtained  to  date  are  presented  in  Tables  1  and  2  for  isotropic 
trapped  protons  and  cosmic  rays  respectively  for  the  geometries  discussed 
above.  The  1-D  case  is  that  studied  in  Refs  2  and  4  but  using  the  revised 
trapped  proton  estimate,  and  in  all  cases  the  nominal  shielding  thickness  of 
0.3  cm  of  ‘LiH  is  used.  It  can  be  seen  that  results  are  in  generally  good 
agreement  between  the  1-D  studies  and  the  3-D  study  on  the  isolated  detector, 
the  latter  case  tending  to  be  higher  by  up  to  a  factor  of  2.  This  is  presumably 
due  to  the  more  accurate  treatment  of  the  cascades  to  include  tertiary  and 
higher  order  products.  Analysis  shows  that  the  production  of  particle  cascades 
enhances  the  spallation  background  over  that  due  to  primaries  by  up  to  a  factor 
of  3  for  trapped  protons  and  up  to  a  factor  of  13  for  cosmic  rays.  The 
additional  production  of  1 2  *1  increases  the  overall  background  enhancement 
factors  to  4.7  and  19.0  respectively.  Comparison  of  the  rates  between  the 
various  geometries  indicates  that  the  presence  of  other  detectors  and  spacecraft 
material  has  a  net  shielding  benefit  for  trapped  protons  but  has  a  deleterious 
effect  for  cosmic  rays.  This  is  due  to  the  greater  penetration  and  higher 
multiplicities  for  the  production  of  secondaries  of  the  more  energetic  cosmic 
rays. 

Using  the  rates  from  the  last  rows  of  each  table  the  best  estimate  of  the 
spectra  of  the  radioactive  background  components  has  been  revised  and  results 
are  presented  in  figures  4  and  5.  The  lower  estimate  for  the  SAA  proton  fluence 
and  the  higher  multiplicities  for  cosmic  ray  effects  raise  the  relative 
importance  of  the  cosmic  ray  components  compared  with  the  estimates  of  Ref  2. 

It  can  be  seen  that  the  production  of  1J,I  by  neutron  capture  provides  a 
significant  component  of  the  background.  Use  of  Monte  Carlo  codes  enables 
neutron  tagging  to  ascertain  the  origin  of  the  neutrons  which  eventually 
capture  in  the  central  Nal  element.  Results  are  presented  in  Table  3  and 
confirm  that  the  majority  of  captures  in  fact  arise  from  neutrons  originating 
in  the  same  detector.  Furthermore,  within  a  detector  the  majority  of  captures 
in  the  central  Nal  are  due  to  neutrons  generated  in  the  surrounding  materials, 
primarily  the  massive  Nal  shield.  This  implies  that  any  neutron  shielding  strategy 
must  aim  to  shield  the  central  Nal  from  the  remainder  of  the  detector  structure. 
However  attempts  to  use  thermalising  absorbers  have  appeared  counterproductive3  . 
This  is  confirmed  in  Figure  6  which  presents  the  influence  of  varying  the 
amount  of  ‘LiH  shielding  upon  12 'I  production  rates  in  the  central  Nal  crystal 
for  the  trapped  proton  case  for  all  the  simulations  and  geometries.  The 
deleterious  influence  of  thermalising  absorber  on  the  l2,I  component  is  confirmed. 
This  is  more  marked  when  there  is  a  large  source  of  neutrons  as  is  the  case  for 
both  the  1-D  study,  when  the  detector  was  allowed  to  grow  with  increased  ‘LiH 
thickness,  and  the  full  spacecraft  simulation.  The  effect  of  an  alternative 
neutron  shield  location  has  also  been  examined.  Naively  it  might  be  supposed 
that  positioning  an  ‘LiH  shield  midway  in  the  outer  Nal  shield  could  provide 
some  advantage  as  it  could  then  thermalise  and  absorb  neutrons  from  the  outer¬ 
most  annulus  of  Nal  while  the  thermalised  neutrons  escaping  would  in  part  be 
removed  by  the  inner  annulus  before  entering  the  central  Nal.  Such  a  design 
would  be  more  likely  to  be  effective  for  the  trapped  proton  source  as  in  this 
ca3e  the  neutron  production  is  more  biased  towards  the  outside.  Hence  a 
simulation  has  been  performed  incorporating  a  2  cm  ‘LiH  layer  at  the  expense 
of  the  central  2  cm  of  the  outer  Nal  sheild.  Results  for  neutron  capture  and 
spallation  are  in  fact  not  significantly  different  (7%  enhanced)  from  the  case 
for  zero  neutron  shield.  Thus  use  of  ‘LiH  shielding  is  not  a  promising  approach 
for  current  sizes  of  detector.  However  it  should  be  noted  that  the  spacecraft 
model  contained  no  thermalising  materials .  If  the  external  neutron  source  has 
a  thermal  component  it  is  possible  that  small  thicknesses  could  have  a  net 
benefit  against  this  component  of  the  background. 
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Preliminary  investigations  have  been  made  of  the  influence  of  detector 
orientation  and  particle  anisotropies.  Reorientation  of  the  detectors  with 
respect  to  the  spacecraft  within  an  isotropic  flux  gives  variations  which  are 
currently  not  statistically  significant  and  at  present  can  be  sized  at  less 
than  15*.  However  use  of  the  trapped  proton  anisotropies  of  Fig  3  as  input 
to  the  simulations  indicates  that  activation  of  the  central  crystal  can  vary 
by  about  a  factor  of  2.3  between  what  are  presumed  to  be  the  best  (spacecraft 
to  the  West  of  the  detectors)  and  worst  (spacecraft  to  the  South  of  the 
detectors,  detectors  pointing  West)  cases  (Fig  7).  Activation  of  the  tungsten 
collimator  is  a  factor  of  4  higher  for  the  worst  case.  Hence  trapped  proton 
effects  could  be  diminished  somewhat  by  appropriate  orientation  during  South 
Atlantic  Anomaly  traversals.  Further  studies  are  needed  using  a  more  accurate 
representation  of  the  collimator  and  covering  a  wider  range  of  orientations 
while  taking  account  of  the  mission  constraints  in  order  to  determine  the  exact 
improvement  that  is  possible. 
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TABLE  1 

TRAPPED  PROTON  EFFECTS 

5  -2 

Daily  trapped  proton  fluence  »  5.0  x  10  cm  ,  Ep >  100  MeV 


Configuration 

|  Mean  Production  Rates  s  '  in  Central  Nal 

Neutrons 

Spallation 

128-j- 

1-D  isolated 

1119 

267 

205 

3-D  isolated 

2033 

552 

266 

3-D  4  detectors  +  10  cm  Al 

1235 

367 

222 

3-D  4  detectors  +  spacecraft 

1038 

378 

248 

TABLE  2 

COSMIC  RAY  EFFECTS 

-2  -1 

Average  Cosmic  ray  flux  *  0.16  cm  s  ,  Ep >  100  MeV 


Configuration 

Mean  Production  Rates  s  1  in  Central  Nal 

Neutrons 

Spallation 

128j. 

1-D  isolated 

1421 

122 

132 

3-D  isolated 

1409 

313 

122 

3-D  4  detectors  +  spacecraft 

1690 

330 

161 

TABLE  3 

ORIGIN  OF  CAPTURED  NEUTRONS 


Proton  Source 

Fraction  of  Total  Captures  in  Single  Detector  by  Neutrons  from: 

1  Same  Detector 

2  Other  Detectors 

3  Spacecraft 

Inner  Belt 

0.73 

0.12 

0.15 

Cosmic  Rays 

0.81 

0.13 

0.06 

Proton  Source 

For  *  1  *  above  Fraction  of  Captures  in  Central  Nal  by  Neutrons  from 

Central  Nal 

|  Csl 

W 

Outer  Nal 

Inner  Belt 

0.16 

0.24 

0.16 

0.44 

Cosmic  Ray 

0.37 

0.13 

0.23 

0.27 

TRAPPED  PROTON  ANISOTROPY  AT  20  MtY 


Figs  1-3 
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Fig  3  Proton  directional  flux  distributions  for  the  SAA  show  a  narrow 
Gaussian  in  cose,  where  e  is  the  pitch  angle,  and  a  preferential  flux 
of  protons  travelling  geomagnetically  eastwards  (  ♦  =  90  ),  particularly 
at  high  energies. 
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Figs  4&5 
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oroic  ac  *ou  *cra»  Results  are  from  the  full  3-D  space-  3-D  simulation  of  spacecraft  plus  detector  is  used  with 

craft  plus  detector  simulation.  isotropic  trapped  proton  fluxes  at  the  revised  intensity. 


Fig  7  Supposed  best  and  worst  case  orientations  with  respect  to  the 
geomagnetic  field  and  the  direction  of  maximum  trapped  proton  fluxes. 
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a  , Carlo  transpor  codes  have  been  employed  to  simulate  particle  transport  of  i osmic 
rays  and  inner)- aelt  protons  in  various  representations  of  the  Gamma  Ray.  Observatory 
Spacecraft  and  he  Orinted  Scintillation  Spectrometer  Experiment.  Results  are  used 
to  accurately  0  .antify  the  contributions  to  the  radioactive  background,  itss'ss 
shielding  options  and  examine  the  effect  of  detector  and  spacecraft  orientation  in 
anisotropic  tr  pped  proton  fluxes. 

V  '  , 


K\(  I  n  .4143  October  I  ‘’SOI 


